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Abstract. Scalar cosmological perturbations of a weakly self-interacting plasma mixed with a 
perfect radiation fluid are investigated. Effects of this plasma are considered through order A 3 / 2 of 
perturbative thermal-field-theory in the radiation dominated universe. The breakdown of thermal 
perturbation theory at vastly subhorizon scales is circumvented by a Pade approximant solution. 
Compared to collisionless plasmas the phase speed and subhorizon damping of the plasma density 
perturbations are changed. An example for a self-interacting thermal field is provided by the 
neutrinos with effective 4-fermion interactions. 

The evolution of cosmological perturbations depends on the matter content of the uni- 
verse. In the radiation dominated epoch photons, electrons, and baryons are coupled strongly, 
hence an effective description as a perfect radiation fluid is possible. The energy density of 
the perfect fluid is related to its pressure by the equation of state p = 3p. The cosmological 
perturbations [[]J are determined by the linearized Einstein equations together with the per- 
turbed equation of state dp = c 2 s 5p, where c s is the sound speed of the perfect radiation fluid; 
the anisotropic pressure vanishes. Neutrinos interact weakly, therefore they decouple from 
the strongly interacting matter at a temperature ~ 1 MeV and propagate freely thereafter. 
This almost collisionless gas of neutrinos may be described by the linearized Einstein- Vlasov 
equations or by a thermal- field-theoretic approach [[J . In the latter the perturbed energy- 
momentum tensor is calculated from the response of the ultrarelativistic plasma on a metric 
perturbation 
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where V is the effective action of (thermal) matter. The Hubble rate H ~ T 2 /mpi anc k 
is much smaller than the temperature T, thus momenta k phys of cosmological interest are 
much smaller than T as well. The second variation of the effective action, the gravitational 
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Figure 1: (a) The co moving density contrast grows as a function of the conformal time r\ on 
superhorizon scales (krj = k phys /H < ir) for the perfect fluid (full line) and the collisionless 
plasma (dotted line). On subhorizon scales (k phys /H > tc) the high sound speed, c s = 1/a/3, 
of the perfect fluid gives rise to acoustic oscillations of constant amplitude. Collisionless 
matter perturbations are damped due to directional dispersion, the phase speed is one, and 
the oscillations are due to the gravitational interactions of the collisionless plasma with itself, 
(b) A mixture of 50% perfect radiation fluid (full line) and 50% collisionless plasma (dotted 
line) evolved together. On subhorizon scales the collisionless plasma feels the gravitational 
potential of the radiation fluid, its damping is modulated. The normalization is arbitrary. 



polarization tensor, is evaluated in the high temperature limit, T 3> fc phys . The thermal-field- 
theoretic formulation goes beyond classical kinetic theory. Eq. ([]]) and the linearized Einstein 
equations define a closed set of integro-differential equations for the metric perturbations, 
which may be solved by a power-series ansatz [[J. Fig. 1(a) shows the comoving density 
contrast 5k as a function of conformal time r\ of a perfect radiation fluid and a collisionless 
plasma respectively. At temperatures below the muon threshold (T < 35 MeV) the ratio a = 
Pu/Ptot ~ 0.49, whereas after e + e~-annihilation at about 0.2 MeV a ~ 0.41. Fig. 1(b) shows 
the density contrasts for a = 1/2. Between temperatures of 1 MeV (decoupling of neutrinos) 
and 10 eV (the cold dark matter contributes only 1/10 of the total energy density) this is 
a good model for the dominant matter perturbations at that time. Neutrinos well above 1 
MeV are included in the radiation fluid usually. How can we describe neutrinos near the 
decoupling temperature and what effects should we expect for cosmological perturbations? 

In a recent work H. Nachbagauer, A. Rebhan, and myself Q studied the behavior of cos- 
mological perturbations for thermal 0(A r )-symmetric scalar fields with weak self -interactions 
A</> 4 . We discovered an unexpected breakdown of thermal perturbation theory on vastly sub- 
horizon scales (k phys 3> H). The full set of 2-loop diagrams, after resummation of thermal 
masses m = A 1//2 T + 0(A) and nonlocal gravitational vertices, has been included in the 
gravitational polarization tensor of Eq. ([!]). The resummation of an infinite subset of all 



2 



0.2 



. 1 




0.05 



0.02 



0.01 



. 005 



. 002 



10 . 



Figure 2: The sub horizon evolution of comoving density perturbations of a collisionless 
plasma (dotted line as in Fig. 1(a)) is compared with the density perturbations of a self- 
interacting scalar plasma with A = 1 (dashed line). At krj/n ~ 4 the thermal perturbation 
theory, evaluated through order A 3 / 2 , starts to break down, which leads to unphysical growth 
of the perturbation asymptotically. The subhorizon behavior of the solutions has been 
improved by use of Pade approximants (full line). Vastly subhorizon the phase speed is 
smaller than one and the damping is weaker than l/(ki]). 

diagrams is standard in hot QCD || and avoids infra-red divergences beyond the leading 
order A. Due to the presence of a thermal mass the next-to- leading order is A 3 / 2 . On su- 
perhorizon scales the behavior of the growing mode is the same for any value of A, but 
the decaying mode may show superhorizon oscillations depending on the values of a and A 
[0, f3j. The subhorizon evolution of the density perturbation is shown in Fig. 2 for A = 1. 
Thermal perturbation theory does not only break down for large values of the coupling A, 
but also vastly subhorizon at perturbatively small values of A. Increasing orders in A 1 / 2 are 
increasingly infra-red singular in the high temperature limit as the external momentum k 
approaches the light-cone. For subhorizon scales rj 3> 1/k the contribution from A 1 overtakes 
the contribution from A , A 3//2 overtakes A 1 , and so on. A way to improve the subhorizon 
behavior is to rewrite the perturbative series in A 1 / 2 into a (2,1) Pade approximant 



In the large N limit this gives an excellent approximation to the exact value of the thermal 
mass |J even for large values of A. A calculation of the most important 1-loop diagram of the 
polarization tensor to all orders in A 1 / 2 has shown that the Pade approximant, in contrast 
to the perturbative result, describes the phases correct and provides a good approximation 
of the damping for large values of A and rj. The full 2-loop Pade improved solution for the 
density contrast with A = 1 is shown in Fig. 2. The main conclusions for scalar perturbations 
(for the vector and tensor sector see ||) are: 

At superhorizon modes the anisotropic pressure decreases as the coupling A is increased. 
At subhorizon scales the phase speed of the damped oscillations is smaller than one, and 
the damping is weaker than l/(kr]). 



a + b\ + cA 3/2 + . . . ~ a + b\/[l - (c/b)X 1/2 ] + ... 
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For one scalar field we have at sub horizon scales the paradoxical situation that, although 
the interaction rate T may be much bigger than the Hubble rate, Y/H ~ A 2 mpi anc k/T, 
when T <C mpi anc k and A is perturbatively small, the plasma perturbations behave almost 
collisionless. This happens because the class of diagrams dropping out in the large N limit 
(order A 2 In A and higher) is not taken into account by the Pade approximant solution. In 
the large N limit, where the Pade approximant is close to the exact solution, the Hubble 
scale picks up a factor \fN from the relativistic degrees of freedom, thus the interactions 
are slower than the expansion. On superhorizon scales the interaction rate is irrelevant for 
the behaviour of the perturbations. The amount of superhorizon anisotropic pressure is 
determined by the initial conditions and gravity only. 

What does this tell us about neutrinos? At low temperatures T 100 GeV the weak 
interactions of neutrinos are well described by effective 4-fermion interactions. The topology 
of the corresponding loop-diagrams is the same as before. Although the details are quite 
different the qualitative features should be similar. Let me naively relate A to GfT 2 . At 
neutrino decoupling this corresponds to A 10~ n , thus perturbation theory should work 
perfect. Our results for the scalar plasma suggest that neutrino perturbations behave almost 
collisionless even above decoupling! On subhorizon scales this is not reliable for two reasons: 
Again our result can only be trusted for a large number of neutrino flavors, and the neutrinos 
couple to matter which is strongly interacting and is described by a perfect fluid. Our results 
show that, depending on initial conditions, effects of anisotropic pressure on superhorizon 
scales can be important although the dominant plasma is collisional. 

I thank H. Nachbagauer and A. Rebhan for their collaboration on the issue of this work. 
The work was supported in part by the FWF and the SNF. 
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